Origin theories for the eccentricities of 
extrasolar planets 



Fathi Namouni^ 

CNRS, Observatoire de la Cote d'Azur, BP 4229, 06304 Nice, France 
namouni@obs-nice . f r 



1 Introduction 

Planetary science has known a revolution in the year 1992 when Wolszczan 
and Frail discovered the first extrasolar planets around a pulsar [69] . This 
discovery was soon followed in 1995 by the first Jupiter-like planet around the 
sun-like star 51 Pegasi [H]. As of this writing we have entered the stage of 
statistics as more than 200 extrasolar planets are known to orbit sun-like stars. 
The revolution that was triggered by these discoveries changed the paradigm 
of what we think of as a planetary system and what we believe to be the 
scenarios that led up to the formation of planets. The planet around 51 Pe- 
gasi is the prototype of what has become to be known as hot Jupiters: planets 
with masses comparable to Jupiter's but with periods of a few days and orbits 
smaller than the tenth of Mercury's distance to the sun. Planetary orbits also 
seem to be eccentric. This observation is not surprising since the solar sys- 
tem planets have eccentric orbits except that half the extrasolar planets have 
eccentricities larger than 0.3 which is much more significant than Jupiter's 
0.05. Such large eccentricities are reminiscent of the small body populations 
of the solar system that got stirred up by their gravitational interactions with 
the larger planets. In this respect, extrasolar planetary eccentricities are more 
unusual, witness the median eccentricity of 0.13 of main belt asteroids larger 
than 50 km. The planetary revolution did not stop at orbital radii and eccen- 
tricities: 10% of known planetary systems belong to binary star systems and 
only one planet so far is in a triple stellar system. Multiple planets around a 
single star make up about 10% of known planets. 

In essence, the planetary revolution heralded the coming of a new plane- 
tary principle: orbital diversity is a rule of planetary formation. Diversity here 
is not meant to imply subtle changes but drastic ones with respect to the as- 
pect of the solar system. This state of affairs has prompted a serious revision 
of the theories of planetary formation: hot Jupiters with few-day orbits could 
not have formed in situ. Instead they have formed outside a few AU for a sun- 
like star where it is possible for ices to condense and for the planets to capture 
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large amount of gas from the protoplanetary disk. Only after they formed, did 
they travel all the way to meet their current orbits. It is interesting to note 
that the concept of radial migration was already known in the contexts of 
accretion disks |32| in binary star systems and of planetary rings |18| . Only 
before 1995, one could not plausibly contemplate the prospect of suggesting 
the existence of massive planets that traveled all the way from Jupiter's cur- 
rent location just to stop on a close orbit with a few-day period. The basic 
aspects of the process of planetary migration through the tidal interaction of a 
planet with the gaseous protoplanetary disk are now well understood [15] yet 
an important challenge remains: what stops planetary migration towards the 
star? The leading contender for stopping planetary migration is the planet's 
interaction with the stellar magnetosphere but a definitive quantitative de- 
scription is still lacking. 

Extrasolar planetary eccentricities have equally resulted in a drastic change 
of perception: it is often heard that it is not the extrasolar planets that are 
eccentric, rather it is the solar system that lacks eccentricity. This perception 
is encouraged by the availability of some simple instabilities that one can set 
up in a many-body gravitational system to simulate the generation of the wild 
orbits of extrasolar planets. Upon close examination such instabilities as well 
as other eccentricity scenarios do not tell the whole story of how extrasolar 
planets become eccentric. In fact, just as the features of the planetary migra- 
tion process yield constraints on the planetary formation scenarios, so do the 
various theories of the eccentricity excitation. 

It is the aim of this chapter to review the various processes of the origin 
of extrasolar planets' eccentricities in the context of planetary formation. We 
start by reviewing the properties of extrasolar planetary orbits in section 2. 
Section 3 contains a commentary on the various known theories of eccentricity 
excitation. Section 4 specializes in a recent addition to the eccentricity theories 
based on an relationship between the planets and the stellar jet that is powered 
by the protoplanetary accretion disk. The final section 5 discusses how the 
eccentricity origin problem may contribute further to the theory of planet 
formation. 



2 Eccentricity observations 

Extrasolar planets are detected with various observational techniques [55j . 
The Doppler analysis of the reflex velocity of the host star is by far the most 
successful technique to date. It is also the technique that has uncovered the 
large eccentricities of extrasolar planets. If a planet has a circular orbit, the 
analyzed stellar spectrum yields a sinusoidal oscillation of the stellar reflex 
motion. If the planet is on an eccentric orbit, the reflex motion as a function 
of time becomes distorted with respect to a pure sine reflecting the unequal 
times the star spends in different locations along its orbit around the center of 
mass of the star-planet system (Figure [1]) . The discovery of large eccentricity 
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orbits by the Doppler reflex velocity method is due to its ability to detect 
planets on wider orbits in contrast to that to the transit method. Planets on 
many-day periods have usually undergone tidal circularization by the host 
star. 
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Fig. 1. Reflex velocity of the stars HD 75289 (left) [63] and HD 108147 (right) 
|51) . The planet around HD 75289 has a circular orbit while that around HD 108147 
has an orbital eccentricity of 0.5. Pictures taken from the Geneva Extrasolar Planet 
Search http://obswww.unige.ch/ ~udry/ planet/ planet.html 



The statistical analysis of extrasolar planet eccentricities reveals very few 
clues as to the origin of the elongated orbits. For the known sample of 196 
planets discovered by the reflex velocity and transit techniques, the median 
eccentricity is at 0.21 if all planets are counted and at 0.28 if planets with 
periods smaller than 5 days are excluded because their circular orbits simply 
reflect tidal circularization. The prevalence of such large eccentricities and 
the large typical mass of the detected planets (comparable to Jupiter's) has 
encouraged the comparison of the extrasolar planetary systems to binary star 
systems. Depending on the methods used, similarities in the eccentricity dis- 
tribution of both populations can be found [5^ or not [21]. What is agreed 
upon is that there is no correlation between the size of the orbits and their 
eccentricities in each population, and no striking resemblance of the scatter of 
both populations in an orbital size versus eccentricity plane. The size of the 
orbit usually refers to either the semi-major axis or the pericentre radius. The 
latter is used to account for those orbits that have not yet had enough time 
to be circularized -as the pericentre distance is conserved under stellar tides. 
Finally, eccentricities show a vague correlation with the planetary masses with 
heavier planets enjoying larger eccentricities. 



4 Fathi Namouni 



3 Eccentricity origin theories 

Seven known explanations have been put forward to account for the large 
eccentricities of extrasolar planets. They are: (1) planet-planet scattering, (2) 
planet-protoplanetary disk interaction, (3) Kozai's secular cycles, (4) excita- 
tion through radial migration into a mean motion resonance, (5) Stellar en- 
counters, (6) stellar-like N-body relaxation, and (7) excitation through stellar 
jet acceleration. In the following, we comment on these possibilities by dis- 
cussing their instability types, characteristic timescales, their epoch of appli- 
cability as well as their advantages and drawbacks. 

3.1 Planet-planet scattering 

Planet-planet scattering is a simple process to generate eccentric orbits in an 
N-body gravitational system. If a system of two or more planets on planar cir- 
cular orbits find themselves "initially" closer than is permitted by Chirikov's 
criterion for the overlapping of mean motion resonances |68| , the planets scat- 
ter off one another leading to a system with more stable albeit eccentric config- 
urations. Depending on the number, masses and "initial" spacings of the plan- 
ets, the instability timescale varies between 10^ to lO'^ years |52 | 167 1 [15 1 140 1 ll6j. 
The epoch that is referred to by the adjective "initial" is that of the disap- 
pearance of the agent or the conditions that kept the planets from scattering 
off one another in the first place. This epoch is customarily associated with a 
significant dispersal of the parent gaseous protoplanetary disk. As well shall 
point out in the next section, planet-disk interaction is known to primarily 
erase orbital eccentricities. An additional condition for planet-planet scatter- 
ing to be operational is the absence of a significant population of smaller 
bodies such as the primordial asteroid belt. Depending on the mass spectrum 
in the planetary system, the smaller populations are able to limit the growth 
of the planetary eccentricities through dynamical friction [4^ . This at least how 
it is believed that the terrestrial planets in the solar system did not acquire 
large eccentricities [U [44] . Numerical works that tackle the extrasolar eccen- 
tricity problem using planet-planet scattering do not consider the effect of 
leftover small-body populations after the gaseous disk has dispersed. Planets 
are set up at a few Hill radii from one another and initial conditions are sam- 
pled to reproduce the eccentricity of certain observed systems. The general 
excitation trend of planet-planet scattering leads to larger eccentricities than 
the ones observed. What may prove to be a serious problem for planet-planet 
scattering is the eccentricity distribution obtained in multiple systems that 
contain Jupiter-mass planets as well as Earth-mass planets. The conservation 
of angular momentum in this case will force the much smaller planets to have 
much larger eccentricities than the Jovian planets. 
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3.2 Planet-disk interaction 

A planet embedded in a gaseous disk excites sound waves at the locations of its 
mean motion resonances within the disk akin to the gravity waves excited by 
Saturn's satellites in its ring system. The density enhancements at the mean 
motion resonances act back on the planet resulting in gravitational torque. 
Two types of resonances contribute to this torque: (1) corotation resonances 
that primarily affect the semi-major axis and tend to damp any acquired ec- 
centricity and (2) Lindblad resonances that primarily affect the eccentricity 
and tend to increase it [T71 [TB] . The torque contribution of the former is larger 
than the latter 's by about 5%. At first sight, planet-disk interaction damps 
the eccentricity on timescale that depends strongly on the disk's thickness and 
less strongly on the disk's mass density and the planet's mass [B11I11I3Z]- The 
torques originating from higher order resonances as well as those pertaining 
to the relative inclination of the planet and the disk do not change the out- 
come significantly [19l |45]. Only if the corotation torque saturates, can the 
Lindblad resonances increase the eccentricity [19l [20] . The conditions under 
which saturation arises are difficult to quantify explaining why an eccentricity 
increase due to a disk-planet interaction has never been observed in numerical 
simulations although this might be due to numerical artefacts [381 . 




3.3 Secular Kozai cycles 

In his study of asteroids perturbed by Jupiter on high eccentricity and incli- 
nation orbits, Kozai [28| showed that the averaging of the interaction poten- 
tial over the mean motion without expanding the force amplitude in terms 
of eccentricity and inclination leads to new types of secular resonances. The 
conservation of the vertical component of angular momentum (vertical refers 
to the direction of Jupiter's orbital normal) shows that when the orbital ec- 
centricity increases, the inclination decreases. In particular, if objects are set 
up on inclined but circular orbits, large eccentricities can be achieved as the 
inclination decreases in its motion around the secular resonance cycle. The 
application of the Kozai cycle to the eccentricities of extrasolar planets as- 
sumes that there is a binary star on a not-too-far inclined orbit that perturbs 
the planet that formed in a circular orbit in a timescale shorter than the 
Kozai libration cycle. In essence, the secular Kozai cycle idea transforms the 
eccentricity problem into an inclination problem. In this sense, the observed 
planets do not possess a proper eccentricity but one that is forced by the stel- 
lar binary and that will always oscillate between its original value, zero, and a 
maximum value depending on the planetary-binary semi-major axis ratio, the 
binary's mass and its orbital inclination with respect to the plane on which 
the planet initially formed. When applied to specific binary star systems with 
one planet, the Kozai mechanism works fine and helps characterize the orbit 
and mass of the secondary star required to excite eccentricity [27l [42l [T4j . Sta- 
tistically, Kozai based excitation of one-planet binary systems tend to yield 
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larger eccentricities than observed [55] ■ As T Tauri stars form in multiple sys- 
tems, it is not unreasonable to try and apply the Kozai mechanism to the 
whole sample of observed extrasolar planets. The problem is that the Kozai 
cycle is usually destroyed by mutual gravitational interactions. The addition 
of more planets to the one-planet binary system, forces the precession of the 
planets' pericentres. If the planets are of comparable mass as it is observed in 
multiplanet systems, the Kozai cycle is lost. 

3.4 Mean-motion resonances 

The role of mean motion resonances in exciting orbital eccentricity has its 
roots in the study of the orbital evolution of Jupiter's and Saturn's regular 
satellites under planetary and satellite tides [50 . These satellite systems are 
known to be in or to have crossed mean motion resonances thereby acquiring 
forced eccentricities. The combined modeling of the orbital evolution, capture 
into resonance and the tidal interaction lends valuable bounds on the dissipa- 
tion factors of Jupiter, Saturn and their satellites. Extrasolar planets form in 
a gaseous disk that does not dissipate after they acquired most of their masses 
or else hot Jupiters would not exist. Planet-disk interaction naturally gives 
rise to orbital migration with different planets in the same system migrating 
at different rates. This differential migration makes planets in the same sys- 
tem encounter mean motion resonances. Capture into resonance may occur 
depending on whether the migration is convergent or divergent (for instance 
if the outer planet is moving faster or slower than the inner planet). Con- 
vergent migration leads to capture into resonance. The subsequent common 
migration of a planetary pair in resonance pumps up the eccentricities on the 
migration timescale [31] [7D] . Divergent migration does not lead to resonance 
capture, instead eccentricity jumps are acquired at resonance passage [7, 62]. 
While convergent migration is certainly the way the known resonant multiple 
systems have acquired their eccentricities, this excitation method involves a 
mystery that may shed light on how to halt planet migration in a gas disk. 
The mystery consists of the observation that convergent migration is far too 
efficient in exciting eccentricities to the point where in many systems, when 
capture occurs, migration must stop quickly thereafter or else eccentricities 
are pumped up to much larger values than those observed. As it is implausi- 
ble to invoke the dispersal of the gas disk, planetary migration may become 
ineffective because of the nonlinear response of the gas disk to the planet 
pair. It is interesting to note that when capture occurs, the planetary relative 
inclinations may be excited as unlike planetary satellites that orbit Jupiter 
and Saturn, the central potential is keplerian. Consequently, for the same or- 
der, eccentricity and inclination resonances are close (but not coincident as 
the gas disk modifies the pericentre and node precession rates). Planet-disk 
interaction is not well understood for large eccentricity planets and off-plane 
(inclined with respect to the disk) orbits. The often used formulas for eccen- 
tricity damping from the disk torques have not been verified for eccentric and 
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inclined planets. Divergent migration has the advantage of being applicable to 
the wider non-resonant multiplanet systems. For a planetary pair, divergent 
migration requires that the inner planet migrates faster than the outer one. 
Gap driven migration (also known as type II) is favorable to such a condition 
as the migration rate is the viscous timescale of the disk. Divergent migration 
may take place because viscosity is likely to be a decreasing function of the 
distance to the star. If however the part of the disk that is located between the 
two planets is dispersed as when the two planetary gaps merge, the direction 
and rate of migration may be altered significantly [39j . 

3.5 Stellar encounters 

Stellar encounters are common events in star clusters. A planetary system that 
encounters a star will have its planets feel a tidal force that elongates their 
orbits . For inner planets that orbit close to the host star, the excitation which 
lasts for about 1000 years will occur on a secular timescale. Outer planets if 
they exist will feel a localized impulse somewhere in their orbits. Typical 
encounter frequencies are of one in 5 x 10^ years while typical encounter 
parameters are a few hundred AU. Unless planets are way outside the classical 
planetary region (inside 30 AU), excitation is not efficient ||71^. To reverse 
this conclusion and account for the eccentricities of inner planets, the system 
must contain several planets with increasing distance and mass from the star 
in order propagate the stellar tug felt by the outermost planet down to the 
innermost ones |71j . 

3.6 Stellar-like relaixation 

The qualitative similarity of the eccentricities of extrasolar planets and stellar 
binaries suggests that planets may form through similar processes as those of 
multiple stellar systems. If planets formed by gravitational instability, the for- 
mation time is so short that the planets find themselves confined to a smaller 
space than their orbital stability permits. The relaxation of such systems leads 
to some planetary ejections and many large eccentricity orbits [49j . The ap- 
plicability of this scenario is limited by two facts: first, the minimum planet 
mass the gravitational instability allows is a few Jupiter masses. This means 
that stellar-like relaxation does not work for planets with masses comparable 
to or smaller than Jupiter's. Second, if a two-phase formation where small 
planets form through rocky core accretion and the larger ones through gravi- 
tational instability [60] , then it is likely that the relaxation of the larger plan- 
ets destroys the smaller planets. This is because the gravitational instability 
timescales are usually smaller than the planetesimal accumulation timescales. 
In fact, if large mass planets form through gravitational instability, they are 
likely to inhibit planetesimal accumulation by clearing the inner disk before 
planetary embryos are born. 
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4 Jet-induced excitation 

Stellar jets enter the eccentricity excitation problem because of their ubiquity 
and simplicity [35]. The story of how this works is as follows: although there 
is disagreement on whether there is a statistically significant resemblance be- 
tween the eccentricity distributions of extrasolar planets and stellar binary 
systems, the qualitative similarity is beyond doubt. Those who wish for the 
similarity to be quantitative, would like to affirm the view that planets are 
the lower end of the outcome of star formation. This question has been settled 
observationally in 2005 with two observations: the first is a hot Saturn with a 
giant rocky core discovered by combining the Doppler reflex velocity method 
with transit photometry |54j . The second is the imaging of the first planetary 
candidate which because of bias due to contrast and resolution happens to 
be a warm distant companion orbiting a young brown dwarf [S]. This proves 
that planets do not need large rocky cores and may form by gravitational 
instability. Exit the link between how planets form and their eccentricities. 

If planets do not form like binary stars, perhaps they undergo similar ex- 
citation processes that lend them similarly elongated orbits. In view of the 
different physical environments where planets and stars form, the simplest 
possible excitation process may depend weakly or not at all on the local dy- 
namics of the stellar or planetary companion. Mathematically, this amounts 
to saying that the acceleration imparted by the process is independent of 
position and velocity. 

Simplicity therefore dictates that the process imparts a constant accelera- 
tion that operates during a finite time window. Simplicity also comes with two 
added advantages: we can already know the excitation time scale and the min- 
imal acceleration amplitude. Dimensional analysis shows that the excitation 
timescale has to be proportional to v/A where v is the keplerian velocity of the 
companion around the main star. Further, if the acceleration is to achieve its 
purpose within the lifetime of the system, v/A must be smaller than about 10^ 
years. This tells us that the acceleration A> 3 x 10^^^('i;/10 kms"^) kms^^. 

The process lacks one more attribute: direction. If the acceleration is inde- 
pendent of the formation processes, its direction cannot depend on anything 
related to the planetary companion such as its orbital plane or the direction 
from the star to the companion. In an inertial frame related to the planetary 
or stellar system, we are not left with much choice but the star's rotation axis. 

To sum up, what we are looking for is a process that appears everywhere 
where planet and star formation takes place, acts like a rocket (i.e. with an 
acceleration that does not depend on the position and velocity of the system) 
and whose direction is related to the star's rotation axis. The answer is stellar 

jets [niig. 

Do planets exist when jets are active? The answer is quite likely. Known hot 
Jupiters have moved close to their host stars because of their interaction with 
the gas disk. So we know the gas disk was present and had viscosity well after 
planets finished forming. The gas accreting on the star because of viscosity 
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is the main ingredient along with the magnetic field that threads it needed 
to launch stellar jets and disk winds. It would therefore be an interesting 
coincidence that jets shut off when planets appear in the gas disk a few AU 
away from the star well outside the jet launching region. 

Are there any observational hints that jet-sustaining disks contain planets? 
The only possible hint so far is the observation of variable brightness asym- 
metries in some jet-sustaining disks ^i56, 6^. The variability timescales of a 
few days to a few years are so small that they imply either a peculiar stellar 
activity in the form of single hot spots or the presence of distortions in the 
disk at the location where the orbital period matches the variablity timescale. 
The first option requires a complex stellar magnetic field that differentiates 
strongly between the two stellar poles. The second option may be caused dis- 
tortions in the disk whose origin could be the presence of embedded compact 
objects. 

Do jets have enough strength to build eccentricity? Inferred mass loss 
rates for known young T Tauri stars lie in the range lO-^Mgyear^i to 
lO~i°M0 year~^ and may be two orders of magnitude larger depending on the 
way the rate is measured from the luminosity of forbidden lines [22l [25l [29] . 
The jet also needs to be asymmetric with respect to the star's equator plane 
or else there would be no acceleration. Interestingly, a growing number of 
bipolar jets from young stars [24l [131 [SOI [35] are known to be asymmetric as 
the velocities of the jet and counterjet differ by about a factor of 2. Mass loss 
processes in young stars therefore yield accelerations: 

A ^ 10-" ( ,] ( A (^] kms-. (1) 

\^10-8Moyear-iy V300kms-V V / 

where M is the stellar mass and Ve is the outflow's high velocity component. 
As jets are time-variable processes, the above estimate is only indicative of 
the epochs at which the rates and velocities are measured. In this sense, it is 
closer to being a lower bound on what accelerations really are over the jet's 
lifetime. What is clear is that asymmetric jet acceleration is larger than the 
minimum amplitude of 10~^^(?;/10 kms~^) kms~^. 

For how long can a jet operate? Jet-induced acceleration is technically no 
different from attaching a rocket to the star and accelerating it very slowly 
with respect to the outer part of the disk and the planets. As a result, the 
star acquires a residual velocity that must be smaller than its orbital velocity 
in the Galaxy or else the star is ejected. In fact, there is an even a stronger 
constraint on the residual velocity from the velocity dispersion in the Galaxy. 
Stars do not have exactly circular orbits in the Galaxy. Their motion is slightly 
distorted or eccentric and such eccentricity is measured as a departure of the 
galactic orbital velocity from that of circular motion. This velocity dispersion 
is known for various stellar populations and is of order a few tens of kilometers 
per second. Since the residual velocity imparted by the jet, V is given by the 
product At where r is the duration of acceleration, imposing that V < (vg) 



10 



Fathi Namouni 



where {vg) is the velocity dispersion in the Galaxy yields: 
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This timescale is shorter than the disk's lifetime. In practice, we shall see that 
shorter times are needed. 

Further excitation properties can be deduced by analyzing the effect of the 
combined jet-induced acceleration and the star's gravitational attraction. As 
the star's pull decreases with distance, there is a specific location where the 
latter matches the jet-induced acceleration (that is independent of position 
and velocity). Outside this radius, the star's pull is weak and orbits escape its 
gravity. This reveals an interesting feature of jet-induced acceleration: stellar 
jets are responsible for the outer truncation of circumstellar disks. It is clear 
that in the interior vicinity of the truncation radius, the orbital perturbations 
are large as the excitation time becomes comparable to the orbital period. In 
this region, the keplerian orbits are subject to a sudden excitation; not only the 
eccentricities are excited but the semi-major axes are also affected leading to 
inward or outward migration. Well inside the truncation radius, the excitation 
time is much larger than the orbital period. In this region, eccentricity builds 
up slowly over a large number of revolutions of the planet around the star 
and the mean orbital radius remains constant on average. Excitation in this 
region occurs on secular timescales. Planetary companions mostly fall inside 
the secular region as they are far inside the truncation radius which is more 
or less the size of the protoplanetary disk. 

4.1 Secular jet-induced excitation 

In the secular region where the excitation time is larger than the orbital period, 
the dynamics of excitation can be simplified by averaging the acceleration 
over the orbital period of the companion. For a constant acceleration, the 
interaction potential is simply R = A ■ x where x is the position vector. 
Averaging the interaction potential amounts to averaging the position vector 
of a pure keplerian motion. A simple calculation shows that (a;) = —3aex{f = 
0)/2r where /, a, e and r are the true anomaly, the semi-major axis, the 
eccentricity and radius of the keplerian orbit. The direction of x at pericentre 
is that of the eccentricity vector e = t; x h/G{m -\- M) — x/\x\ where v is 
the velocity vector of the companion, G is the gravitational constant, m and 
M are the masses of the companion and the host star and h = x x v is the 
specific angular momentum. This enables us to write the secular potential as: 



where in the last equality, the ^-direction of the reference frame is chosen along 
along A and w, f2, /, are the longitude of pericentre, longitude of ascending 
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node and the inclination of the orbit. To simphfy the excitation problem 
further, we use the conservation of the component of angular momentum 
h along the direction of acceleration as A ■ h = A ■ {x x A) = 0. In the 
reference frame where A is along the z-direction, the conservation of angular 
momentum yields (1 — e^)^/^ cos / = cos Iq where Iq is the initial inclination of 
the keplerian orbit with respect to the jet-induced acceleration. This relation 
enables us to eliminate the inclination variable in {R) and reduce the problem 
to an integrable, one dimensional system with: 

{R) = --Aad ^^^^ esm^. (4) 

where w = ti7 — i7 is the argument of pericentre. The time evolution of eccen- 
tricity and argument of pericentre is obtained from: 

d{R) . d{R) 
^ = 2 a — ' 2 

where n = ^J G{M + m) /a^ is the companion's mean motion. In this one- 
dimensional system, e and uj follow curves of constant (i?) shown in Figure 
([2]). There are equilibria at uj — ±90° and e = •\/2sin(/o/2) corresponding 
to 7 = cos~^(v'cos/o). The maximum value of e is sin/o and corresponds 
to the cycle of initially circular orbits. For these orbits, (i?) = throughout 
their cycle implying that the orbits orientation can take only one value uj — 
modulo 180°. 
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Fig. 2. Contour plots of the acceleration potential (|4|) in the eccentricity e and 
argument of pericentre uj{°) plane. The direction of acceleration makes an angle 
lo ~ 30° with respect to the companion's angular momentum vector. The time 
evolution of the two orbits (e = 0, a; = 0) and (e = 0.3, ui = 90°) is shown in Figure 
©■ 
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For time-dependent accelerations and provided that the variation timescale 
is longer than the orbital period, the eccentricity evolution is given by: 

. 3A{t)e 



2na 



sin^/o — e^, (6) 



where e is the sign of cos a; which is set by the requirement that e > 0. The 
solution of ^ can be found exactly as: 
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sin /o 



(7) 



The inclination is obtained from cos / = cos /q / — e(r)2. For strictly con- 
stant accelerations (infinite time window), A{t) = Aq and e oscillates between 
and sin Jo at the excitation frequency: 

3|^o| 

UA = . (8) 

na 

Examples of such oscillations that were obtained from the direct integration 
of the full equations of motion are shown in Figure ([3]). The good agreement 
between the secular solution and the results of the numerical integration comes 
from the fact that A is independent of position and velocity. 

To optimize the excitation of a finite eccentricity from an initially circular 
state, the duration of acceleration needs to be smaller than half the oscillation 
period: t < 7rna/3|Ao|. Examples of eccentricity excitation at three different 
semi-major axes (i.e. three different excitation frequencies) are shown in Fig- 
ured The dependence of the excitation amplitude on the ratio of the duration 
to the excitation time is also seen in the same Figure. Note that because ec- 
centricity excitation in the secular region is a slow process compared to the 
orbital time, the convolution of the dynamics under the constant acceleration 
A with a finite time window has the effect of shutting off the excitation at 
some eccentricity value depending on the duration. 

The secular excitation through jet-induced acceleration is therefore able 
to make e reach sin Iq and is largest if the initial orbital plane contains the 
direction of acceleration {Iq = 90°). As (i?) = for initially circular orbits, cu 
and fi, remain at zero. This forcing of the pericentre to be perpendicular to 
the direction of acceleration favors apsidal alignment in multiplanet system. 

If the jet's inclination with respect to the companion's orbital plane is small 
Jq <C 1, the maximum eccentricity will be negligible. The sin/o limitation is 
problematic because it is not reasonable to expect stellar jets to be highly 
inclined with respect to the gas disk where the companions form. Fortunately 
for the jet-induced excitation theory, there is a natural way out. Jets are 
known to precess over timescales from 10^ to 10"* years [lOl [HI [59] . The origin 
of such precession is not known as we lack resolution to probe inside the jet 
launching region. It is possible that precession is caused by a warp in the 
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Fig. 3. Time evolution under a conservative acceleration. The eccentricity e, ar- 
gument of pericentre and inclination /(°) are shown for an initially circular 
orbit e = (solid) and an orbit librating about the secular resonance uj = 90° with 
an initial eccentricity e = 0.3 (dashed). The semi-major axis is identical for both 
orbits and is set to unity. The acceleration corresponds to a period of 10* years at 
1 AU. The plots were obtained by the numerical integration of the full equations of 
motion. 



disk's plane resulting from interactions with stellar companions as T Tauri 
stars are known to form in multiple systems. Precession is attractive because 
it offers the possibility of resonance if the excitation frequency ua matches the 
jet precession frequency J7a- This in fact is exactly what happens when the 
eccentricity evolution is derived in the situation where the constant magnitude 
acceleration rotates at a constant rate. It turns out that the corresponding 
secular problem is also integrable. The eccentricity and inclination evolution 
are given by [IS] : 

2 p^sin^a ,- 9\ ,/H X ,/, X 

e = 5 — ^ 2(3 +p ) — 4(1 +pcosa) cosj/+i — 4(1 — pcosa) cosi/_i 

+ (1 - + V+V-) cos{v+ - v^)t +(1 - - ly+iy-) cos{i'+ + z/_)i] (9) 
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Fig. 4. Eccentricity excitation by time-dependent constant-direction accelerations. 
The equations of motion are integrated numerically with an acceleration A{t) — 
AoH{t) exp —t/r where H{t) is the Heaviside unit step function and Ao = 2.21 x 
10~^^kms~^. The oscillation period at 1 AU is 1.11 x 10^ years. The timescale 
r = 7200 is chosen so that V =5kms^^. The curves correspond the semi-major 
axes: 1 AU (sohd), 32 AU (dashed) and 128 AU (dotted). 

cos/ = J- [(/ - p2 + 2 +p2[p2 _ 3] cos2a) 

zi/^v^Vl — 

+p'^ sin (p^ + 1 + 2p cos a) cos i^+t (10) 
-l-p^ sin (p^ + 1 - 2pcosQ;) cosi^_t] (11) 

where a is the jet angle with respect to the z-axis of the reference frame, 
Vj- = + 1^2p cos a, p = ua/^SIa, and the time t is normalized by S7a- The 
companion's initial orbit is circular and lies in the xy-plane. 

Nominal resonance is defined where the frequency match, p — I, occurs. 
It corresponds to a nominal resonant semi-major axis aios given as: 

(^) ( 2xlO-.'l.....-0 "(Tfei)" ™ 

where Tpicc = Itx j flA- In terms of the resonant semi-major axis, the frequency 
ratio can be written as p = \/ a/ arcs- Far inside resonance {p ^ 1), the jet 
precesses faster than the eccentricity excitation leading to a reduction of the 
eccentricity amplitude from sin a to 2psina. Far outside resonance (p S> 1), 
the jet's precession is slow compared to the eccentricity excitation so that 
the latter is described by a constant acceleration without rotation. In the 
resonance region, the proximity of p to unity increases the denominators of 
the eccentricity expression ([SJ which leads to eccentricities close to unity. At 
exact resonance, the eccentricity reaches unity regardless of the jet angle. The 
width of the region around resonance increases with the jet angle a. These 
features are illustrated in Figures (O where we plot the expressions (jH]) and 
(fTTj) for a jet angle a ~ 1°, an excitation time I-k juA = 10'* years, and the 
three values of p: 0.05, 0.9, and 1. Finally, we note that as the eccentricity 
excitation time is ua-, no resonant forcing occurs when Qa = nvn the secular 
region (jia ^ n). 
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Fig. 5. Excitation of the eccentricity (solid) and inclination (dashed) by a nearly 
perpendicular precessing jet with an angle a = 1°. The companion's orbit is located 
at a = 1 AU and evolves from a circular orbit in a plane orthogonal to the jet's 
precession axis. The acceleration is Ao = 2 x 10~^''kms~^ yielding an excitation 
time of 2-K/nA = 10* years. From top to bottom, the panels were obtained from 
equations (|9l jlip with the frequency ratios, p: 0.05, 0.9, and 1 - the precession period 
is 2p X 10* years. The symbols correspond to the numerical integration of the full 
equations of motion. 



4.2 Sudden jet-induced excitation and radial migration 

The location where orbits can not longer be retained by the star is where 
the frequency ua becomes comparable to the local mean motion n of the 
companion. Near this limit, the forced periodic oscillations of the semi-major 
axis a are reinforced by the eccentricity and acquire large amplitudes making 
the orbits unstable in the long term. Calling Okpir the semi-major axis of 
the keplerian boundary where the star where ua — n, we the jet-induced 
acceleration is expressed as: 
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corresponding to an excitation period Ta = 27r/n^: 

Figure ([5]) shows an example of an escape orbit of a constant-direction accel- 
eration with akpir = 10^ AU and an inclination Iq — 30°. The orbit's initial 
semi-major axis is 68.5 AU. The characteristics of motion are not strictly ke- 
plerian as the companion hovers above the star. Such escape orbits offer an 
interesting way to expel planets from around their parent stars or equivalently 
to disrupt a binary stellar system. If a companion is formed near the keple- 
rian boundary or is pushed out to it by a possibly remaining inner disk that 
followed photo-evaporation |26| . it could become unbound. 
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Fig. 6. Escape of a companion located near the keplerian boundary of a conservative 
constant-direction acceleration with Jo = 30°. The distances are given in AU. The 
boundary's semi-major axis is flkpir = 100 AU. Note how the companion hovers 
above the star before escaping. 

For a realistic jet, the induced acceleration has a finite duration. Accord- 
ingly, akpir varies in time from infinity before the jet's launch to a location 
determined by the strongest acceleration the jet can provide. Ultimately, the 
keplerian boundary is pushed out to infinity. The acceleration's finite dura- 
tion extends the keplerian boundary depending on the ratio of the duration r 
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to the excitation time at the keplerian boundary of the equivalent constant- 
acceleration problem T4(akpir)/2. When r > r4(akpir)/2, orbits beyond Okpir 
have enough time to acquire sufficient momentum and escape the gravita- 
tional pull of the star. When r < TA(akpii)/2, the stability region extends 
beyond a^pir- The new stability boundary is given by the semi-major axis 
where r ~ TA(aoo)/2 which is larger than Okpir since Ta is a decreasing func- 
tion of the semi-major axis a. The expressions of Ta and the residual velocity 
V, show that Uoc — G{M + m)V^'^, the location where the keplerian velocity 
V matches V. 

Orbits near the keplerian boundary of a finite duration acceleration that 
do not escape the pull of the star will end up with elongated orbits whose 
semi-major axes and eccentricity have changed. This happens because the 
companion feels an almost instantaneous velocity kick (the orbital period is 
large compared to r) . The conservation of linear momentum and energy can 
be combined to find the change in semi-major axis as: 

1 _ 1 2Fsin/oCose' ^2 

o^^o^ '^/G(mT^^, G{M + m) ' 

where Oi, Of are the initial and final semi-major axis, 9 is the longitude of the 
companion along its orbit, and /q is the inclination of the orbital plane with 
respect to the direction of the residual velocity V. Note that for /q ^ 0, the 
final semi-major axis can be larger or smaller than the initial value depending 
on the longitude of the companion where the velocity pulse if felt. 




1e+02 1e+03 ai 1e+04 1e+02 1e+03 ai 1e+04 



Fig. 7. Migration and eccentricity excitation near the keplerian boundary of a finite 
duration acceleration for two inclination values 7o = (left panel) and Jo = 20° 
(right panel) . In each panel, the final semi- major axis at ( AU) (solid) and the final 
eccentricity a (dashed) are shown as a function of the initial semi-major axis ai (AU). 
The parameters are: V — 0.35 km s~^, akpir = 300 AU, and r — 500 years. The 
inclined circular orbits were started at the descending node {6 — 180°). For Jo = 0°, 
aoo = 7341 AU and for 7o = 20°, = 14 542 AU. 



Figure ([7]) shows the final semi-major axes and eccentricities at two differ- 
ent inclinations Iq = and 20° for an acceleration amplitude corresponding to 
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Okpir = 300 AU and a duration r — 500 years resulting in a residual velocity 
y = 0.35 km s^^. The companions were started at the descending node for 
Jq = 20° to illustrate inward and outward migration. Such migration could 
enhance the delivery of minor bodies to the Oort Cloud and explain the trans- 
port of Kuiper Belt outliers 2000 CR105 and Sedna (90377). 

4.3 A test case: the v Andromedae binary system 

Multiplanet systems provide good test cases for the excitation by jet-induced 
acceleration. This is because mutual gravitational interactions cause the ec- 
centricity vectors to precess. If the ensuing precession rates are much faster 
than the excitation frequencies, the slow build up of eccentricity by jet-induced 
acceleration will be lost [46( 47, • This situation is similar to that encountered 
for precessing jets inside the resonance radius Ores- 

An interesting system for testing the excitation mechanism is that of v An- 
dromedae [3J. It contains three planets, two of which have their apsidal direc- 
tions aligned [53l|33l[6] as well as a 0.2Mo stellar companion at a projected 
distance of 750 AU [36]. 

Numerical simulation can be used to reproduce the planetary orbits from 
initially circular co-planar orbits to their observed state [13]: af, — 0.059 AU, 
ef, = 0.020, = 241°, mbsini = 0.75 mj, = 0.821 AU, = 0.185, = 
214°, mcsini = 2.25 mj, aa = 2.57 AU, = 0.269, = 247°, and mdsini = 
2.57 mj. It turns out that mutual planetary perturbations are strong enough to 
prevent excitation if the acceleration is smaller than Aq ~ 10~^^kms~^. The 
equivalent smallest keplerian boundary is at akpir ^ 500 AU. Below this value, 
the current configuration can be recovered along with the apsidal alignment 
of the outer two planets. The presence of the stellar companion outside the 
keplerian boundary leaves us two options: either the excitation by acceleration 
is ruled out or that the companion was initially inside the boundary and 
migrated during a sudden excitation (the projected distance of 750 AU does 
not translate necessarily into a semi-major axis as the companion's orbit is 
likely to be eccentric). 

Figure [8] shows a simulation of the jet-induced acceleration of the form 
A(t) = Aq/ cosh{t ~ to)/T where 3 x 10~"kms"^ and r = 2000 years, 

applied the current planetary system plus a stellar companion on an orbit 
of semi-major axis a = 298 AU, just inside the keplerian boundary of A, 
cikpir ~ 300 AU. The stellar companion's initial orbit was given an eccentric- 
ity e = 0.3 in order to decouple its perturbations from the planetary sys- 
tem. In particular, the eccentricity excitation by the Kozai mechanism [281 is 
not efficient because the corresponding excitation time (~ 10^ years) is much 
larger the duration of acceleration and the eccentricity secular frequency of 
the isolated two-planet system (~ 7000 years). The jet-induced acceleration 
produces a configuration similar to the observed one with stellar orbital el- 
ements: e — 0.5, a = 600 AU. Apsidal alignment is achieved as the result of 
the acceleration's strength that maintain the forcing of companion orbits to 
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Fig. 8. Eccentricity excitation, apsidal alignment and binary migration in the v And 
system. The acceleration pulse is shown not to scale in all panels. The first panel 
shows the eccentricity excitation of planets d (solid) and c (dashed) and the ec- 
centricity evolution of v And B (short-dashed) as well as the acceleration pulse 
normalized to its maximum value (dotted). The second panel shows the radial mi- 
gration of V And B. The last panel shows the relative apsidal libration of planets d 
and c. 

be perpendicular to the direction of acceleration. Note that only when the 
acceleration's strength is near maximum and the keplerian boundary nears 
300 AU, does the stellar orbit acquire a larger eccentricity. 

4.4 The solar system 

Was the solar system subject to jet-induced acceleration? There are two ob- 
servations that hint at the dynamical action of the solar system's jet. The first 
is the inclination of Jupiter's orbital normal by 6 degrees with respect to the 
sun's rotation axis. As Jupiter is the more massive planet in the solar system, 
this implies that either the early protoplanetary disk of the solar system was 
warped with respect to the Sun's equator plane or that the Jupiter's orbit 
gained inclination with respect to an early equatorial disk. Both possibilities 
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are consistent with the jet-induced acceleration theory. More recently, the dis- 
covery of Calcium- Aluminum inclusions in the grains of comet Wild 2 suggest 
that the only possibility for the Jupiter-family comet originating from the 
early Kuiper belt to contain such high temperature minerals is that they were 
transported by the solar system's jet. An additional piece of the puzzle is given 
by outer solar system bodies such as Sedna (90377) that are decoupled from 
the solar system's planets as their perihelia are larger than Neptune's orbital 
radius. Such objects could have been transported by the radial migration in 
the sudden excitation region associated with the jet-induced acceleration. 

Why are the solar system's planets eccentricities small? There are three 
possible reasons for this: first, the resonance radius ^2]) where the jet's preces- 
sion matches the excitation time could have been outside the planets' orbits. 
Inside the resonance radius, very little excitation can take place (Figure [5]). 
Second, mutual planetary perturbations could have destroyed the secular ec- 
centricity growth. Third, if the jet angle were to be small and the acceleration 
duration equally brief, the planets' location near the resonant radius would 
be of little help to raise their eccentricities. 

4.5 The unknowns of jet-induced excitation 

The unknowns of jet-induced excitation belong to two categories: (1) the un- 
knowns of disk-planct interaction and (2) the unknowns of the time evolution 
of jets. In the first category comes the issue of the eccentricity damping by 
the accretion disk. As explained in the section about mean motion resonances, 
the damping of eccentricity for large eccentricity and inclination orbits is not 
currently understood. The advantage of jet-induced excitation with respect 
to the excitation during migration while in mean motion resonance is that 
excitation times can be much smaller than the migration timescale and the 
viscous timescale of the disk. In this phase, the formed planets gain a substan- 
tial inclination that makes them exit the gas disk which should in principle 
reduce the eccentricity damping significantly. In this context, it is useful to 
bear in mind that jet-induced acceleration becomes effective in a planetesimal 
disk only when a few planets are left. The planetesimal mutual gravitational 
interactions destroy excitation through the radom precession of their orbits. 
When the bodies left in the disk are such that the precession periods due to 
their mutual perturbations are larger than the excitation time, jet-induced 
acceleration becomes effective and planets may exit the disk on inclined or- 
bits. For the second category of unknowns, it is safe to say that except for 
precession, little information is available about the time variations of jets over 
their entire life span. Perhaps the main advantage of jet-induced excitation is 
its small set of parameters: amplitude, duration and jet precession frequency. 
These determine all the features of eccentricity growth or lack thereof. In 
multiplanet systems, the acceleration subjects all companions, planetary and 
stellar alike, to the same instability as it is independent of position and veloc- 
ity. A better knowledge of the time dependence of acceleration can therefore 
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easily confirm or rule out the effect of jet-induced excitation in multiplanet 
systems. 

5 Concluding remarks 

The various theories of eccentricity excitation are valuable tools to gain insight 
into one of the most pressing if not the most pressing issue in planet formation 
theory: the mismatch of protoplanetary disk lifetimes and the timescales of 
planet formation and planet migration. Protoplanetary disk lifetimes range 
from 10^ to 10^ years. To have a viable theory, the planet formation timescale 
has to be smaller by at least an order of magnitude than these estimates. 
Planets are believed to form in two ways: (1) rocky core accretion and (2) 
gravitational instability. The formation of rocky cores precedes gas accumu- 
lation on the way to forming gaseous planets. This specific phase is slow. 
Typical formation timescales range from 10® to 10^ years at Jupiter's location 
|34| . Gravitational instability however has a much smaller timescale of order 
10'^ years at the same location. With all its caveats about the required disk 
opacity for gravitational fragmentation and the role of shearing instabilities 
in disrupting a forming protoplanet, formation by gravitational instability 
looked far more promising than rocky core accretion to explain the existence 
of the hot Jupiters. This was true until the discovery of a hot Saturn with a 
70 Earth mass rocky core orbiting the star HD 149026 [54] -for reference, the 
mass of Jupiter's rocky core is believed to be 15 Earth masses. The forma- 
tion of a 70 Earth mass rocky core either as a single body or as a merger of 
smaller cores is likely to last at least 10^ years according to current planetes- 
imal accumulation models. This temporal crisis does not stop at rocky core 
formation, it worsens because of planet migration. Migration arises from the 
planet's interaction with the gas disk and has typical timescales of 10® years 
for an Earth mass core (linear regime, type I) to 10^ years for a Jupiter mass 
planet (nonlinear regime, type II) 65 , 9J. The duration of the eccentricity ex- 
citation phase and its dependence on the mechanism's parameters may help 
elucidate the time sequence of the events that produced the observed extra- 
solar planets. In particular, three-dimensional hydrodynamic simulations of 
disk-planet interactions for eccentric and inclined orbits (off the disk's mid- 
plane) may elucidate the problem of the end of resonant excitation in multiple 
systems. Planet-planet scattering models may benefit from including a wider 
mass spectrum in the populations of simulated planets. The excitation by 
Kozai's secular mechanism requires further assessments of the ability to form 
multiple planets under the stellar companion's perturbations [61^ . Similarly, 
jet-induced excitation was demonstrated for a single planet under the action 
of a precessing acceleration. The effect of mutual planetary entrainment for 
planets on both sides of the precession-excitation resonance needs to be in- 
vestigated and applied to observed extrasolar systems. 
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